Tomography of silicate dust around M-type AGB stars I. Diagnostics based on dynamical models

The heavy mass loss observed in evolved asymptotic giant branch stars is usually attributed to a two-step process: atmospheric levitation by pulsation-induced shock waves, followed by radiative acceleration of newly formed dust grains. Detailed wind models suggest that the outflows of M-type AGB stars may be triggered by photon scattering on Fe-free silicates with grain sizes of about 0.1 - 1 $\mu$m. Due to the low grain temperature, these Fe-free silicates can condense close to the star, but they do not produce the characteristic mid-IR features that are often observed in M-type AGB stars. However, it is probable that the silicate grains are gradually enriched with Fe as they move away from the star, to a degree where the grain temperature stays below the sublimation temperature, but is high enough to produce emission features. We investigate whether differences in grain temperature in the inner wind region, which are related to changes in the grain composition, can be detected with current interferometric techniques, in order to put constraints on the wind mechanism. To investigate this we use radial structures of the atmosphere and wind of an M-type AGB star, produced with the 1D radiation-hydrodynamical code DARWIN. The spectral energy distribution is found to be a poor indicator of different temperature profiles and therefore is not a good tool for distinguishing different scenarios of changing grain composition. However, spatially resolved interferometric observations have promising potential. They show signatures even for Fe-free silicates (found at 2-3 stellar radii), in contrast to the spectral energy distribution. Observations with baselines that probe spatial scales of about 4 stellar radii and beyond are suitable for tracing changes in grain composition, since this is where effects of Fe enrichment should be found.Comment: Accepted for publication in Section 8. Stellar atmospheres of Astronomy and Astrophysics. The official date of acceptance is 07/09/2017. 9 pages, 7 figures, 4 figures in appendi

Exploring wind-driving dust species in cool luminous giants III. Wind models for M-type AGB stars: dynamic and photometric properties

Stellar winds observed in asymptotic giant branch (AGB) stars are usually attributed to a combination of stellar pulsations and radiation pressure on dust. Shock waves triggered by pulsations propagate through the atmosphere, compressing the gas and lifting it to cooler regions, which create favourable conditions for grain growth. If sufficient radiative acceleration is exerted on the newly formed grains through absorption or scattering of stellar photons, an outflow can be triggered. Strong candidates for wind-driving dust species in M-type AGB stars are magnesium silicates (Mg$_2$SiO$_4$ and MgSiO$_3$). Such grains can form close to the stellar surface, they consist of abundant materials and, if they grow to sizes comparable to the wavelength of the stellar flux maximum, they experience strong acceleration by photon scattering. We use a frequency-dependent radiation-hydrodynamics code with a detailed description for the growth of Mg$_2$SiO$_4$ grains to calculate the first extensive set of time-dependent wind models for M-type AGB stars. The resulting wind properties, visual and near-IR photometry and mid-IR spectra are compared with observations.We show that the models can produce outflows for a wide range of stellar parameters. We also demonstrate that they reproduce observed mass-loss rates and wind velocities, as well as visual and near-IR photometry. However, the current models do not show the characteristic silicate features at 10 and 18 $\mu$m as a result of the cool temperature of Mg$_2$SiO$_4$ grains in the wind. Including a small amount of Fe in the grains further out in the circumstellar envelope will increase the grain temperature and result in pronounced silicate features, without significantly affecting the photometry in the visual and near-IR wavelength regions.Comment: 11 pages, 14 figure

Molecular opacities for low-mass metal-poor AGB stars undergoing the Third Dredge Up

The concomitant overabundances of C, N and s-process elements are commonly ascribed to the complex interplay of nucleosynthesis, mixing and mass loss taking place in Asymptotic Giant Branch stars. At low metallicity, the enhancement of C and/or N may be up to 1000 times larger than the original iron content and significantly affects the stellar structure and its evolution. For this reason, the interpretation of the already available and still growing amount of data concerning C-rich metal-poor stars belonging to our Galaxy as well as to dwarf spheroidal galaxies would require reliable AGB stellar models for low and very low metallicities. In this paper we address the question of calculation and use of appropriate opacity coefficients, which take into account the C enhancement caused by the third dredge up. A possible N enhancement, caused by the cool bottom process or by the engulfment of protons into the convective zone generated by a thermal pulse and the subsequent huge third dredge up, is also considered. Basing on up-to-date stellar models, we illustrate the changes induced by the use of these opacity on the physical and chemical properties expected for these stars.Comment: 23 pages, 8 figures, accepted for publication in Ap

An extensive grid of DARWIN models for M-type AGB stars I. Mass-loss rates and other properties of dust-driven winds

The purpose of this work is to present an extensive grid of dynamical atmosphere and wind models for M-type AGB stars, covering a wide range of relevant stellar parameters. We used the DARWIN code, which includes frequency-dependent radiation-hydrodynamics and a time-dependent description of dust condensation and evaporation, to simulate the dynamical atmosphere. The wind-driving mechanism is photon scattering on submicron-sized Mg$_2$SiO$_4$ grains. The grid consists of $\sim4000$ models, with luminosities from $L_\star=890\,{\mathrm{L}}_\odot$ to $L_\star=40000\,{\mathrm{L}}_\odot$ and effective temperatures from 2200K to 3400K. For the first time different current stellar masses are explored with M-type DARWIN models, ranging from 0.75M$_\odot$ to 3M$_\odot$. The modelling results are radial atmospheric structures, dynamical properties such as mass-loss rates and wind velocities, and dust properties (e.g. grain sizes, dust-to-gas ratios, and degree of condensed Si). We find that the mass-loss rates of the models correlate strongly with luminosity. They also correlate with the ratio $L_*/M_*$: increasing $L_*/M_*$ by an order of magnitude increases the mass-loss rates by about three orders of magnitude, which may naturally create a superwind regime in evolution models. There is, however, no discernible trend of mass-loss rate with effective temperature, in contrast to what is found for C-type AGB stars. We also find that the mass-loss rates level off at luminosities higher than $\sim14000\,{\mathrm{L}}_\odot$, and consequently at pulsation periods longer than $\sim800$ days. The final grain radii range from 0.25 micron to 0.6 micron. The amount of condensed Si is typically between 10% and 40%, with gas-to-dust mass ratios between 500 and 4000.Comment: Accepted to A&A, 17 pages, 15 figure

Exploring wind-driving dust species in cool luminous giants II. Constraints from photometry of M-type AGB stars

The heavy mass loss observed in evolved asymptotic giant branch (AGB) stars is usually attributed to a two-stage process: atmospheric levitation by pulsation-induced shock waves, followed by radiative acceleration of newly formed dust grains. The dust transfers momentum to the surrounding gas through collisions and thereby triggers a general outflow. Radiation-hydrodynamical models of M-type AGB stars suggest that these winds can be driven by photon scattering -- in contrast to absorption -- on Fe-free silicate grains of sizes 0.1--1\,$\mu$m. In this paper we study photometric constraints for wind-driving dust species in M-type AGB stars, as part of an ongoing effort to identify likely candidates among the grain materials observed in circumstellar envelopes. To investigate the scenario of stellar winds driven by photon scattering on dust, and to explore how different optical and chemical properties of wind-driving dust species affect photometry we focus on two sets of dynamical models atmospheres: (i) models using a detailed description for the growth of Mg$_2$SiO$_4$ grains, taking into account both scattering and absorption cross-sections when calculating the radiative acceleration, and (ii) models using a parameterized dust description, constructed to represent different chemical and optical dust properties. By comparing synthetic photometry from these two sets of models to observations of M-type AGB stars we can provide constraints on the properties of wind-driving dust species. Photometry from wind models with a detailed description for the growth of Mg$_2$SiO$_4$ grains reproduces well both the values and the time-dependent behavior of observations of M-type AGB stars, providing further support for the scenario of winds driven by photon scattering on dust.Comment: Accepted for publication in A&A. 15 pages, 14 figure

Modelling the atmosphere of the carbon-rich Mira RU Vir

Context. We study the atmosphere of the carbon-rich Mira RU Vir using the mid-infrared high spatial resolution interferometric observations from VLTI/MIDI. Aims. The aim of this work is to analyse the atmosphere of the carbon-rich Mira RU Vir, with state of the art models, in this way deepening the knowledge of the dynamic processes at work in carbon-rich Miras. Methods. We compare spectro-photometric and interferometric measurements of this carbon-rich Mira AGB star, with the predictions of different kinds of modelling approaches (hydrostatic model atmospheres plus MOD-More Of Dusty, self-consistent dynamic model atmospheres). A geometric model fitting tool is used for a first interpretation of the interferometric data. Results. The results show that a joint use of different kind of observations (photometry, spectroscopy, interferometry) is essential to shed light on the structure of the atmosphere of a carbon-rich Mira. The dynamic model atmospheres fit well the ISO spectrum in the wavelength range {\lambda} = [2.9, 25.0] {\mu}m. Nevertheless, a discrepancy is noticeable both in the SED (visible), and in the visibilities (shape and level). A possible explanation are intra-/inter-cycle variations in the dynamic model atmospheres as well as in the observations. The presence of a companion star and/or a disk or a decrease of mass loss within the last few hundred years cannot be excluded but are considered unlikely.Comment: 15 pages. Accepted in A&

The complex environment of the bright carbon star TX Psc as probed by spectro-astrometry

Context: Stars on the asymptotic giant branch (AGB) show broad evidence of inhomogeneous atmospheres and circumstellar envelopes. These have been studied by a variety of methods on various angular scales. In this paper we explore the envelope of the well-studied carbon star TX Psc by the technique of spectro-astrometry. Aims: We explore the potential of this method for detecting asymmetries around AGB stars. Methods:We obtained CRIRES observations of several CO $\Delta$v=1 lines near 4.6 $\mu$m and HCN lines near 3 $\mu$m in 2010 and 2013. These were then searched for spectro-astrometric signatures. For the interpretation of the results, we used simple simulated observations. Results: Several lines show significant photocentre shifts with a clear dependence on position angle. In all cases, tilde-shaped signatures are found where the positive and negative shifts (at PA 0deg) are associated with blue and weaker red components of the lines. The shifts can be modelled with a bright blob 70 mas to 210 mas south of the star with a flux of several percent of the photospheric flux. We estimate a lower limit of the blob temperature of 1000 K. The blob may be related to a mass ejection as found for AGB stars or red supergiants. We also consider the scenario of a companion object. Conclusions: Although there is clear spectro-astrometric evidence of a rather prominent structure near TX Psc, it does not seem to relate to the other evidence of asymmetries, so no definite explanation can be given. Our data thus underline the very complex structure of the environment of this star, but further observations that sample the angular scales out to a few hundred milli-arcseconds are needed to get a clearer picture

Catching the fish - Constraining stellar parameters for TX Psc using spectro-interferometric observations

Stellar parameter determination is a challenging task when dealing with galactic giant stars. The combination of different investigation techniques has proven to be a promising approach. We analyse archive spectra obtained with the Short-Wavelength-Spectrometer (SWS) onboard of ISO, and new interferometric observations from the Very Large Telescope MID-infrared Interferometric instrument (VLTI/MIDI) of a very well studied carbon-rich giant: TX Psc. The aim of this work is to determine stellar parameters using spectroscopy and interferometry. The observations are used to constrain the model atmosphere, and eventually the stellar evolutionary model in the region where the tracks map the beginning of the carbon star sequence. Two different approaches are used to determine stellar parameters: (i) the 'classic' interferometric approach where the effective temperature is fixed by using the angular diameter in the N-band (from interferometry) and the apparent bolometric magnitude; (ii) parameters are obtained by fitting a grid of state-of-the-art hydrostatic models to spectroscopic and interferometric observations. We find a good agreement between the parameters of the two methods. The effective temperature and luminosity clearly place TX Psc in the carbon-rich AGB star domain in the H-R-diagram. Current evolutionary tracks suggest that TX Psc became a C-star just recently, which means that the star is still in a 'quiet' phase compared to the subsequent strong-wind regime. This is in agreement with the C/O ratio being only slightly larger than 1.Comment: 11 pages, 9 figures, 5 table

Abundance analysis for long period variables. Velocity effects studied with O-rich dynamic model atmospheres

(abbreviated) Measuring the surface abundances of AGB stars is an important tool for studying the effects of nucleosynthesis and mixing in the interior of low- to intermediate mass stars during their final evolutionary phases. The atmospheres of AGB stars can be strongly affected by stellar pulsation and the development of a stellar wind, though, and the abundance determination of these objects should therefore be based on dynamic model atmospheres. We investigate the effects of stellar pulsation and mass loss on the appearance of selected spectral features (line profiles, line intensities) and on the derived elemental abundances by performing a systematic comparison of hydrostatic and dynamic model atmospheres. High-resolution synthetic spectra in the near infrared range were calculated based on two dynamic model atmospheres (at various phases during the pulsation cycle) as well as a grid of hydrostatic COMARCS models. Equivalent widths of a selection of atomic and molecular lines were derived in both cases and compared with each other. In the case of the dynamic models, the equivalent widths of all investigated features vary over the pulsation cycle. A consistent reproduction of the derived variations with a set of hydrostatic models is not possible, but several individual phases and spectral features can be reproduced well with the help of specific hydrostatic atmospheric models. In addition, we show that the variations in equivalent width that we found on the basis of the adopted dynamic model atmospheres agree qualitatively with observational results for the Mira R Cas over its light cycle. The findings of our modelling form a starting point to deal with the problem of abundance determination in strongly dynamic AGB stars (i.e., long-period variables).Comment: 13 pages, 22 figures, accepted for publication in A&